Superconductivity-Induced Anderson Localisation. 
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We bave studied tbe efFect of a random superconducting 
order parameter on tbe locabzation of quasi-particles, by nu- 
merical finite size scabng of tbe Bogobubov-de Gennes tigbt- 
binding Hamibonian. Anderson locabzation is obtained in 
d = 2 and a mobility edge wbere tbe states localize is observed 
in d = 3. Tbe critical bebavior and localization exponent are 
universal witbin error bars botb for real and complex random 
order parameter. Experimentally tbese results imply a sup- 
pression of tbe electronic contribution to tbermal transport 
from states above tbe bulk energy gap. 
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During the past few years phase-coherent transport in 
hybrid superconducting structures bas emerged as a new 
field of study, bringing together the hitherto separate ar- 
eas of superconductivity and mesoscopic physics. Récent 
experiments have revealed a variety of unexpected phe- 
nomena, including zero-bias anomahes ||ï|-^, re-entrant 
and long-range behaviour Q and phase-periodic trans- 
port l^-^]. Thèse experiments can ail be described 
by combining traditional quasi-classical Green's function 
techniques with boundary conditions derived initially by 
Zaitsev pÔ[ | and simplified by Kuprianov and Lukichev 
[ pT[ or alternatively by generalised current-voltage re- 
lations Q based on a multiple scattering approach to 
phase-coherent transport. The latter approach focusses 
attention to Andreev scattering p3| , whereby an électron 
can coherently evolve into a hole and vice versa, without 
phase breaking. 

The aim if this Letter is to address a new phenomenon, 
not describable by quasi-classical techniques, namely the 
onset of quasi-particle Anderson localisation due to spa- 
tial fluctuations in a superconducting order parameter. 
In contrast with ail of the above experiments, wbere the 
superconducting order parameter A(r) is typically homo- 
geneous, there are many situations in which A(r) varies 
randomly in space, even though the undcrlying normal 
potential is perfectly ordered. One example is provided 
by the melting of a flux lattice [|l4|jl5[ in an otherwise 
perfectly crystalline high Te superconductor. Another 
should occur in anisotropic superconductors, where by 
analogy with ^He-A, disordered textures can arise when 
an anisotropic phase is nucleated from a more symmet- 
ric phase such as '^He-B. In the first of thèse examples, 
the order parameter is not quenched. Nevertheless, close 



to the melting curve, the time scale for changes in A(r) 
can be made arbitrarily long and therefore in the spirit 
of the Born - Oppenheimer approximation, it is reason- 
able to freeze the disorder and when necessary, treat any 
temporal fluctuations as a contribution to the inelastic 
scattering lifetime. 

In one dimension, it is straightforward to demonstrate 
Ip^ that fluctuations in A(r) alone can localise the exci- 
tations, even at énergies high above the bulk energy gap. 
However, localisation in strictly one-dimension is of little 
interest experimentally and therefore in this Letter, we 
address the question of whether or not superconductivity 
induced Anderson localisation occurs in higher dimen- 
sions. Early analytic work by |ï^,|ï^, suggested that in 
the présence of time reversai symmetry, states of energy 
E — are localised for dimensions d < 2, while in the ab- 
sence of time reversai symmetry such states are localised 
in ail dimensions. However calculations using a numeri- 
cal finite size scaling approach ||ï^ were inconclusive and 
to date there has been no expérimental confirmation of 
thèse prédictions. In this Letter we provide the first firm 
numerical évidence for superconductivity induced Ander- 
son localisation in d = 2 and 3 dimensions and for the 
first time compute the exponent v controUing the diver- 
gence of the localisation length ^ at the mobility edge in 
d= 3. 

To address the question of superconductivity induced 
Anderson localisation, we analyze the tight-binding 
Bogoliubov-de Gennes équations 
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Ec^,{E) = ^eME) + 1* ^j{E) + A*ME) 

where ipi{E) {(j)i{E)) indicates the particle (hole) wave- 
function of energy E on site i and j sums over the neigh- 
bours of i. Since only scaling behaviour near a criti- 
cal point is of interest, we examine the simplest possible 
model of a System with no normal disorder, but a spa- 
tially fluctuating order parameter, obtained by choosing 
Éi equal to a constant eo for ail sites i and to set the 
energy scale, choose j — 1. Two models of disorder 
will be examined. In model 1, (which préserves time 
reversai symmetry), we choose A^ — Ao[l + S Ai] and 
in model 2, (which breaks time reversai symmetry), we 
choose Aj = Ao[(l + S Ai) + + (5A-], where S Ai and 
ôA'^ are random numbers uniformly distributed between 
—SA and +SA. In what foUows, we choose eo = 0. 
For each model, we compute the transfer matrix T for 
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a long strip {d = 2) and a long bar {d = 3) of length 
L sites and cross-section M''"^ sites, respectively and 
identify the localisation length with the inverse of 
the corresponding smallest Lyapunov exponent. The re- 
sults are, of course, sensitive to the chosen energy E and 
since, in the absence of disorder (ie 6 A = 0), there ex- 
ists an energy-gap at i? = 0, the usual choice of i? = 
adopted in the absence of superconductivity is inappro- 
priate. As a guide to a reasonable choice of E, we con- 
sider the related problem of a System with normal disor- 
der but with a uniform order parameter. In this case is 
chosen randomly from a uniform probability distribution 
but A; = Aq, for ail i. As noted in ||§ if ip'^iEo) is a 
solution of 



E=A. e.=0 8a=0.1 
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FIG. 1. (a) The /M plotted as a function of the finite width 
M in two dimensions for E = Aq, eq = 0, SA = 0.1, and varions 
values of Aq. (b) As in (a) but in three dimensions where a critical 
point is indicated. 



the normal-state Schrodinger equa- 

tion, namely e,^P°iEo) ~ V'j(^o) = Eo,p°{Eo), then 
ipi{E) and (j)i{E) are each proportional to ijj^{E()), where 



E = \/[Eq + |Aop]. This means that if in the absence 
of superconductivity a state at energy Eq is localised by 
normal disorder, then in the présence of a uniform or- 
der parameter Aq, quasi-particle states at energy E are 
localised with the same localisation length. As a con- 
séquence ail critical properties are unchanged, provided 
Eq is replaced by E. In the normal state problem the 
least localised states occur at Eq = and therefore in 
the présence of normal disorder and a uniform supercon- 
ducting order parameter thèse states occur at E = |Ao|. 
Of course, in what foUows we are interested in the oppo- 
site limit of a spatially fluctuating order parameter with 
no normal disorder. Nevertheless, guided by the above 
observation we choose E =< |Ai| >, where < |Ai| > 
is the ensemble averaged order parameter, which gives 
E — Aq for model 1 and E = v^Aq for model 2. 

The raw data for ^m/M versus M, for model 1 with 
-E = Ao, eo = and SA = 0.1, arc shown in figures 
l(a) and (b) for two and three dimensions, respectively. 
The strength of disorder in the order parameter is VF = 
2Ao(5A, whose critical value is denoted Wc, and is varied 
by changing Aq, with fixed SA. In two dimensions /M 
decreases with increasing M indicating that ail states are 
localised with Wc — 0, whereas in three dimensions there 
is a cross-over from localised to extended behaviour at 
around Aq « 12 which for the adopted value of SA — 0.1 
corresponds to Wc « 2.4. 

To quantify the critical behaviour in 3 dimensions, we 
lincarizc the data about Wc by writing \og{£^M /M) — 
c^m+Pm log W and obtain the coefficients «m and (3m for 
varions M. In terms of the fixed point values logi^M /M)c 
and log Wc, we note that a^/ = log(ÇM /M)c — (3m log Wc- 
Thus, a graph of a^/ versus Pm yields \og{£,M /M)c, 
— log Wc and hence the critical disorder Wc- The critical 
exponent v for the divergence of the localization length 
^ of the infinité System is obtained by substituting aM 
into the first linear relation, which yields log{^M /AI) — 
log(ÇM/M)c + PM^og{W/Wc)- Moreover, near the crit- 
ical point log{W/Wc) (W - Wc)/Wc and Ç - jW - 
Wc\-\ so that log(6//M) = logi^M/M)c ± T'/'/Îm, 
where the -I- (-) sign refers to W > Wc (W < Wc)- The 
finite size scaling requirement S,m/M = f{^/M) immedi- 
ately implies Pm ~ M^^'^, which permits the computa- 
tion of the exponent v- 

Figure 2 shows a graph of log{S^M /M) versus logAg, 
from which a m and /3m for the chosen widths M can 
be extracted. The top-right insert shows the result- 
ing plot of aM versus Pm whose slope is — log Wc and 
the corresponding intercept is log(^jv//-^^)c- This yields 
Wc = 2.36 ±0.04 which corresponds to Aq 11.73±0.12 
and {^m/M)c = 0.58 ± 0.02. The lower-left insert shows 
log /3m versus log M whose slope yields the critical expo- 
nent = 1.64 ± 0.06. 



2 



E=A, 5A=0.1 




guish such small différence if it exists. 
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FIG. 2. Log-log plot of Çm/A"/ versus Aq where the intersec- 
tion defines Wc and {Çm /M)c- The upper-right insert shows the 
coefficients a m versus fSj^j whose slope is — log Wc ■ The lower-left 
insert shows a log-log plot of /3jv/ versus M which yields the value 
for the exponent u = 1.64 it 0.06. 



For model 2, where time reversai invariance is broken 
due to the présence of a complex order parameter, ail 
States are localized in d —2. In contrast, figure 3(a) shows 
the corresponding plots of £,m/M versus M in 3 dimen- 
sions which clearly show a cross-over from extended to 
localised behaviour. Results from a more accurate cal- 
culation are presented in figure 3(b), where the upper- 
right figure yields Wc = 5.57 ± 0.12 and {^m/M)^ = 
0.58 ± 0.02. The lower-left insert shows log/3jif versus 
logAf, the slope of which leads to the value for the ex- 
ponent 1^ = 1.69 ± 0.06. The errors in the calculation 
of ^m/M are monitored as a function of length L and 
chosen to be less than about 0.01 by taking long strips of 
lengths L = 250000 and bars of L = 200000 (L = 50000) 
for the real (complex) case. The errors for Wc and v 
are estimated from the corresponding least-square fits. 
We have also repeated our calculations by taking points 
doser to the critical value W^ where the above analysis 
holds, with no significant change of our results. 

The first important feature of the above calculation is 
the unambigous prédiction of superconductivity-induced 
quasi-particle localization in d =2 and the présence of 
a mobility edge in d —i. Localization arises from fluc- 
tuations in the superconducting order parameter alone, 
without the need for additional normal disorder. A sec- 
ond key resuit is the observation that for both models 
we find {£,m/M)c ~ 0.58 and v ^ 1.6, which are re- 
markably close to the values reported for normal d =3 
real Systems [Q, and also consistent with reported data 
for ordinary disordered critical Systems with and with- 
out time-reversal invariance p^ , |23t . Recently, sligthly 
différent scaling behavior is obtained with and without 
time-reversal by an alternative data analysis based on 
polynomial fits p4|. Our study in d = 3 cannot distin- 
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FIG. 3. (a) Log-log plot of Çm/A'/ versus M for varions values of 
Aq when time reversai symmetry is broken (model 2) which shows 
a cross-over from extended to localised states. (b) Log-log plot of 
S;m/M versus Aq for model 2 where the intersections define Wc 
and {£,M /M)c. The upper-right insert shows the coefficinets anf 
versus /Bm whose slope is — log Wc- The lower-left insert shows a 
log-log plot of /3jv/ versus M which yields the value for the exponent 
u = 1.69 ± 0.06. 



From an expérimental point of view, it is worth noting 
that the absence of quasi-particle diffusion does not imply 
the vanishing of the electrical conductance, because An- 
dreev scattering does not conserve quasi-particle charge. 
It does, however, imply a vanishing of the electronic con- 
tribution to thermal transport from certain states above 
the gap. In a clean superconductor at a finite tempéra- 
ture T, this varies as exp(— A/fcf,T), where A is the bulk 
energy gap. In contrast, in the présence of a fluctuation- 
induced quasi-particle mobility edge Ec, this will be re- 
placed by exjp(—Ec/khT). Thus, for example, the melt- 
ing of a flux lattice in a high température superconductor 
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should be accompanied by an cxponcntial change in the 
electronic contibution to the thermal conductance. 
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